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Abstract 

The method of QCD sum rules in the presence of the external electromagnetic 
field is used to analyze radiative decays of charmed or bottomed mesons 
such as D* — > and B* — > B^y, with the susceptibilities obtained previously 
from the study of baryon magnetic moments. Our predictions on D* decays 
agree very well with the experimental data. There are differences among the 
various theoretical predictions on B* decays but the data are not yet available. 



PACS Indices: 14.40.Jz; 13.20.Fc; 13.20. Jf 



Heavy flavor physics provides a natural and convenient laboratory for testing basic 
ideas such as the CKM matrix, heavy quark symmetry, and aspects related to quantum 
chromodynamics (QCD). Indeed, it is often believed that the nonperturbative aspects of 
QCD may play a less significant role in the heavy quark sector as compared to that in 
the light quark systems, thereby offering simplifications in testing fundamental ideas. 

The situation regarding radiative decays of charmed and bottomed mesons such as 
D* — ► D7 and B* — ► B'-f is a little confusing at the moment. In the charmed sector, 
the predictions [1,2] based on heavy quark symmetry differ from those [3] obtained by 
implementing simultaneously both heavy quark symmetry and chiral symmetry. And, 
yet, these predictions also differ significantly from those f| obtained in a relativistic 
constituent quark model. The disparity among these predictions makes it very difficult, 
if not impossible, to assess the significance of the underlying idea(s) under investigation. 

As a potentially useful benchmark, we wish to present and discuss the predictions 
which may be obtained by employing the method of QCD sum rules [[J, especially using 
such method in the presence of an external electromagnetic field || 0]. The important 
aspect here is that the various susceptibilities involved in such method have already been 
determined from previous studies of baryon magnetic moments and as a result our pre- 
dictions are parameter-free on top of the straightforward nature of applying the method 
to the heavy quark sector. 

We recall || |7| that the two-point correlation function H(p) in the presence of the 
constant electromagnetic field may be used as the basis for the derivation of QCD 
sum rules. Using the radiative decay B*° — > B ^ as a specific example, we have 

n„(p) = i / d 4 xe^(0\T[jS* °(x), j B °(O)]\O) Fa0 , (1) 

where the interpolating fields are 

J'f V) = b^x^d^x) , j B \x) = T{x)i lb d a {x) , (2) 
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with a the color index. The decay constants Jb* and fs are defined as 

(0|jf (0)|S*(p)) = / B . m B » e M , 



(0|j*(0)|B(p)) = / 



(3) 



' m b +m d ' 

where is the polarization 4-vector of the .B* meson. 

The decay process £>*(1~) — > -B(0~) + 7 may be characterized uniquely as follows: 

£ = e?^e tll/XS &'A»B* x d s B , (4) 

T7lB* 

where is the electromagnetic field. The corresponding partial width is 

T(B*^B 1 ) = a^^^, (5) 
3 m B * 

with Pcm the i? meson momentum in the CM frame (i.e. in the rest frame of the initial 
vector meson B*). 

The external field F^ u may induce changes in the physical vacuum and modify the 
propagation of quarks. Up to dimension six (d < 6), we may introduce three induced 
condensates as follows. 

(0\qcr^q\0) F ^ = e q exF^(0\qq\0) , 
9c(0\q^G; u q\0) F ^ = e q eKF^(0\qq\0) , (6) 
g c e^(0\q^G^q\0) F ^ = ie q e^F^ (0\qq\0) , 

where q refers mainly to u or d with suitable modifications for the s quark (while the 
condensates involving the charm or bottom quark are assumed to be negligible), e is the 
unit charge, e u = |, and = — ~. In studying baryon magnetic moments, Ioffe and 
Smilga H obtained the nucleon anomalous magnetic moments \i v = 3.0 and fi n = —2.0 
(±10%) with the quark condensate susceptibilities \ ~ ~~ 8 GeV 2 and k, £ quite small. 
Balitsky and Yung |7j adopted the one-pole approximation to estimate the susceptibilities 
and obtained, 

x = -3.3 GeV -2 , k = 0.22, f = -0.44 . (7) 
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Subsequently Belyaev and Kogan || extended the calculation and obtained an improved 
estimate x — —5.7 GeV -2 using the two-pole approximation. Chiu et al [ID also estimated 
the susceptibilities with two-pole model and obtained 

X = -4.4 GeV" 2 , « = 0.4, £ = -0.8. (8) 

Furthermore, Chiu and co-workers (JToJ re-analysed the various sum rules by treating Xi 
k, and £ as free parameters to provide an overall fit to the observed baryon magnetic 
moments. The opitimal values which they obtained are given by 

x = -3GeV" 2 , k = 0.75, £ = -1.5. (9) 

We note that in these analyses the susceptibility values are consistent with one another 
except for the earliest result x = —8 GeV -2 in ||, which is considerably larger (in mag- 
nitude). In what follows, we shall adopt the condensate parameters x = — 3.5GeV~ 2 , 
2k + £ pd with k, £ quite small, which represent the average of the latter values dis- 
cussed above. The same set of susceptibilities yields a very good overall agreement || 
with the experiment for the vector meson decays in SUf(3) sector (it, d, s). 

At the quark level, the correlation function n^fj)) may be evaluated with the aid of 
the formula: 

n» = -%J d 4 xe ipx Tr[tS e b f (-x)^iS f d e (x)i l5 } (10) 

where iS ab (x) is the quark propagator in the presence of the external eletro magnetic field 
F^v 0. To incorporate the (heavy) quark mass, we choose to work in momentum space. 
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For example, we have, for the down quark propagator, 
iSfiv) = 5 ab ^- 

a \r J p—^d 

+\ h f9cGl vW ^ w [a^ip + m d ) + (p + m d )a^)} 

+ le d 5 ab F^ w ^[a^(p + m d ) + (p + m d )a^)\ 

-5 ab ^(2n) 4 5 4 (p) 

-5 ab ^^(27r) 4 g^d^dJ 4 (p) 

-S ab e d ^[a ■ Fg^ - \fa ■ F-f]^)%d v S^) 

-S ab e dX ^a ■ F(2n) 4 S 4 (p) 

-5 ab e d K^[a ■ F(T ~ \fa ■ FY](2n)%dJ 4 (p) 



+ ■■•, 

with p = jo M 7 M . Here we follow |6|, [7]] and do not introduce induced condensates of higher 
dimensions. 

Substituting the quark propagator back into n M (p) (and neglecting condensates re- 
lated to the heavy quark), we obtain 

+ t r d (o|^|o>F^ + ¥(*-«- i)(o|^Io>f4f (12) 

+ f|(0|^a • GdlO)^^ - (§e 6 + f )(0|^ • Gd|0> 
+ • • -K/^i^V 

Note that the operators with lowest dimension (d = 1) arise from the heavy quark mass 
rrih and the induced light quark condensate x(0|gg|0), which are respectively the first two 



terms in Eq. (|12|) . The next two terms represent dimension-three contributions arising 
from the three condensates (0|gg|0), K(0\qq\0), and C(0|gg|0) . At dimension five, we have 
kept only the quark-gluon mixed condensate [the remaining two terms in ([12]) ] as the 
contributions from gluon condensates are found to be negligible. 

The short- distance operator-product expansion (|I2]) converges rapidly. For bottomed 
mesons, the perturbative term contributes up to as much as 90% of the total sum rules 
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(at p 2 = 0), while the induced quark condensate dominates the power corrections and its 
contribution is about 9%. The contributions from the condensates of higher dimensions 
cancel among themselves and the net results are very small. For charmed mesons, the 
perturbative contribution is about 70% as due to the comparatively smaller charm quark 
mass. Again the contribution from the induced quark condensate is dominant among the 
remaning power corrections. We wish to point out that, unlike in the QCD sum rules 
analyses of charmonium where the perturbative contribution and the gluon condensate 
play a very important role, in the present analysis the contributions from the light quark 
condensates x(0|gg|0) and (0|gg|0) are in fact much more important than that from the 
gluon condensate. 

At the hadronic level, we express H(p) as follows: 

n^-(p) = / b ./b m } 9B : B " 2 1 2 2 1 5 e^F*V + • • ■ , (13) 
p rrib + ma 2 p z — m B , p z — m B 

where we write down explicitly the leading B* — > B'-f contribution and denote the contri- 
bution from the excited states and continuum simply by the ellipsis. 

Unlike for light quark systems, Borel transform is no longer the most efficient way to 
accentuate the contribution from the ground state. The reason is twofold: First of all, we 
note that in performing Borel transform the single-pole contributions from excited states 
give rise to an additional term AM B on the phenomenological side of the sum rule [H, TA . 



Due to the small energy gap between the ground state and excited states in heavy meson 
systems, this term could be fairly large. On the other hand, it was suggested |TT] that 



a method of taking suitable moments (as used below) could be more efficient to suppress 
the contributions from the excited states, as the heavy quark mass provides a natural 
scale (m h > A QCD ). 

To simplify the numerical manipulation we adopt the following approximation: 

1 1^1 

p 2 — m 2 B * p 2 — m 2 B (p 2 — m 2 B ) 2 ' 

where m B = m B*+ m B _ approximation is justified to the level of one part in 10 4 (since 



m>B* — rriB = 46 MeV and we are very far away from either of the two poles). 

In order to pick out the ground state contribution we define moments by taking 
derivatives of the correlator H(p) (the coefficient of e^^F^p 17 in U^p)) P, |TT{], 

Mn = ^(^rn(p)\ p2= _ Q ., (Qg >o). (is) 

At the phenomenological side (r.h.s.) we have 



where 5 n {Qo) contains the contributions from the higher resonances and the continuum. 
For high n, 5 n (Ql) goes to zero rapidly because of the factors + ft? 2 ) n+2 it contains and 
because % > %. So from a certain n onwards M£ h ' s '(Qo) will be practically equal to 
the contribution of the ground state. We also introduce the ratios of the moments: 

(0 2 ) = M n ^(Q 2 ) _ 2 2 i + <y w _i(gg) 

which immediately yields the mass if we are at sufficiently high n where S n (Qo) ~ 
^n-i(Qo) [0- We use the ratios to check the convergency and consistency of our sum 
rule. As observed in |5|, 11 1 the stability region in n will change with the parameter 



o 

y = and one should study the stability of the moments for growing y. A plateau 
develops for certain n and y if the sum rules work well, which arises from the balance 
of the pertuabative contributions and power corrections. We have checked that, in the 



region y = 0.3 ~ 0.7 and n = 11 ~ 16, the mass % = y — ^ — - extracted in this way 
agrees very well with the experimental value. The heavy quark mass is scale dependent. 
We take rrij, = 4.7 GeV and m c = 1.3 GeV as commonly adopted in the QCD sum rule 
approach WM. We take the values from the standard two-point sum rules for the decay 



constants and /#* [14 



f B = 140 MeV , fs* = 160 MeV . (18) 



Using the experimental heavy meson masses mg = 5.28 GeV and mg. = 5.33 GeV [12 
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and equating 

M l n h - s iQl)=M r n h - s iQl) (19) 

with y, n in the region y — 0.3 ~ 0.7, n — 11 ~ 16 where the plateau develops, we finally 
obtain 

g B «o BOy = -4.0 ±0.4, 

g B * +B+1 = 6.8 ± 0.6 , (20) 
9b*b s1 = -5.0 ± 0.5 . 

We have included the correction due to the strange quark mass m s = 150 MeV which 
reproduces roughly the B* — B* mass splitting. Note that the indicated errors are due 
to the method of QCD sum rules itself. Typically in the region y — 0.3 ~ 0.7 and 
n — 11 ~ 16 the perturbative term, the induced quark condensate, the dimension-three 
correction, and the dimension-five correction contribute 56%, 57%, —18%, and 5% of the 
sum rule, respectively. Using the above couplings, we obtain the partial width as follows: 

r B .°-B<> 7 = 0.13 ± 0.03keV , 

rV + ^B +7 = 0.38±0.06keV, (21) 
T B *^ Bsl = 0.22 ± 0.04keV , 

which are also listed in Table I. The experiments to measure these partial widths are yet 
to be carried out. 

The above derivation holds for charmed mesons, except that the charm quark is not 
so heavy as the bottom quark so that the plateau develops at larger y (y = 1.5 ~ 2.0, 
n = 5 ~ 8). Using the sum rule values f D = 170 ± 10 MeV and f D , = 240 ± 20 MeV and 
the observed masses m D = 1.87 GeV and ran* = 2.01GeV, we obtain the couplings: 

fi , _D*°D°7 = 2.88 ± 0.3 , 
<7d*+d +7 = -0-38 ±0.1, (22) 
= -0.3 ±0.1. 

The sum rule for g D *o DOry is stable. However, the sum rule for g D *+ D + 1 suffers from the 
fact that the perturbative contributions from the charm and up quarks tend to cancel 
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each other. To extract gn*+D+-y re hablely, we combine the sum rules for g D *+ D + 1 and 
fi f D*°D°7 to get a new stable sum rules for the quantity (gD*+D+-y + 9D*°D°-y)- Nevertheless, 
the error is larger for go*+D+-y The partial widths as calculated from the above couplings 
are given below. 

r D „o^ D o 7 = 12.9 ± 2keV , 
r D »+^ D+7 = 0.23±0.1keV, (23) 
r D .^ D „ = 0.13±0.05keV, 
which are also listed in Table I. 



Note that there are data from the CLEO collaboration flT5|: 



0.572 ±0.057 ±0.081, 



S = = 0.04 ± 0.06 ± 0.06 , (24) 

1 i D*°^I)0 7 



B(D* + -> D+ 7 ) = 0.011 ± 0.014 ± 0.016 . 



Our results yield 



0.018, (25) 

7 

which is consistent with the CLEO data. On the other hand, the partial width for 



D* — > D ± 7T has also been obtained making use of the light-cone QCD sum rules [14 . 

T(D*° -> D°n°) = (23 ± 3)keV , (26) 

T(D* + -> D°7T + ) ± r(£>* + -> D+7T ) = (46 ± 7)keV . (27) 
Combining our predictions with these results, we obtain 

rVo^ D o 7 _ q 5g (2g) 



1 D ,o^ £) o 7r o 

£(£>*+ -> D+ 7 ) = 0.005 , (29) 
both of which are in agreement with the CLEO data cited above. 
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TABLE I. Comparison between some theoretical calculations based on heavy quark and 
chiral symmetry 0,0], the relativistic quark model ||, and QCD sum rules approach (this 
work). The unit is keV. 





Cho a § 


Cho 6 i 


Cheng § 


Jaus [g] 


this work 


D*° - 


-> D° 7 


11.3 ±22.0 


8.8 ± 17.1 


34 


21.69 


12.9 ±2 


D* + - 


-> D+7 


10.6 ± 10.5 


8.3 ±8.1 


2 


0.56 


0.23 ±0.1 


D\ - 








0.3 




0.13 ±0.05 


B*° - 


-> S°7 


0.163 ±0.261 


0.127 ±0.203 


0.28 


0.142 


0.13 ±0.03 


B*+ - 


-> S+7 


0.846 ± 1.494 


0.66 ±0.931 


0.84 


0.429 


0.38 ±0.06 


B* s - 


^7 










0.22 ±0.04 



a with the parameters m c = 1500 MeV, = 4500 MeV. 
b with the parameters m c = 1700 MeV, = 5000 MeV. 

To sum up, we collect in Table I the existing theoretical predictions on the radiative 
decay widths for the bottomed and charmed mesons. As indicated earlier, our predictions 
in the charm sector are in good overall agreement with the existing CLEO data. In the 
bottom sector, there are differences (albeit less significant) among the various theoretical 
predictions but the data are not yet available. 
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